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[i]  We  report  the  first  evidence  of  artificial  ionospheric 
plasmas  reaching  sufficient  density  to  sustain  interaction 
with  a  high-power  HF  pump  beam  produced  by  the 
3.6  MW  High-Frequency  Active  Auroral  Program  (HAARP) 
transmitter  in  Gakona,  Alaska.  The  HF-driven  ionization 
process  is  initiated  near  the  2nd  electron  gyroharmonic 
at  220  km  altitude  in  the  ionospheric  F  region.  Once  the 
artificial  plasma  reaches  sufficient  density  to  support 
interaction  with  the  transmitter  beam  it  rapidly  descends  as 
an  ionization  wave  to  —150  km  altitude.  Although  these 
initial  artificial  layers  appear  to  be  dynamic  and  highly 
structured,  this  new  ability  to  produce  significant  artificial 
plasma  in  the  upper  atmosphere  opens  the  door  to  a  new 
regime  in  ionospheric  radio  wave  propagation  where 
transmitter-produced  plasmas  dominate  over  the  natural 
ionospheric  plasma  and  may  eventually  be  employed  as 
active  components  of  communications,  radar,  and  other 
systems.  Citation;  Pedersen,  T.,  B.  Gustavsson,  E.  Mishin, 
E.  Kendall,  T.  Mills,  H.  C.  Carlson,  and  A.  L.  Snyder 
(2010),  Creation  of  artificial  ionospheric  layers  using  high-power 
HF  waves,  Geophys.  Res .  Lett.,  37 ,  L02106,  d o i ;  10. 1 029/ 
2009GL041 895. 

1.  Introduction 

[2]  The  ionospheric  plasma  a  few  hundred  kilometers 
above  the  earth’s  surface  acts  as  a  reflector  for  radio  waves 
in  the  HF  frequency  range,  enabling  long-distance  radio 
communications,  over-the-horizon  radar,  and  other  technol¬ 
ogies  of  growing  importance  for  national  security  and  counter¬ 
terrorism  applications  [ McNamara ,  1991].  However,  the 
natural  diurnal,  seasonal,  and  solar-driven  variability  of  the 
ionosphere  imposes  severe  limitations  on  operation  of  such 
systems.  Concepts  to  actively  control  radio  propagation  by 
using  high-power  electromagnetic  waves  to  ionize  the  gas¬ 
ses  in  the  upper  atmosphere  and  create  artificial  ionospheric 
plasma  have  been  discussed  for  decades  [e.g.,  Koert ,  1991], 
but  have  never  been  demonstrated  due  to  the  large  electric 
field  thresholds  required  for  conventional  or  even  runaway 
breakdown  [Gurevich  and  Zybin ,  2005].  Recent  advances 
in  high-power  HF  excitation  of  the  ionosphere  made  possi¬ 
ble  by  completion  of  the  High-Frequency  Active  Auroral 
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Research  Program  (HAARP)  transmitter  facility,  however, 
have  made  production  of  significant  artificial  ionospheric 
plasma  possible  for  the  first  time,  opening  the  door  to  prac¬ 
tical  exploration  of  these  concepts. 

[3]  Electromagnetic  waves  with  ordinary  (O)  mode  polar¬ 
ization  are  reflected  from  plasmas  at  the  point  the  trans¬ 
mitted  wave  frequency  fj  equals  the  local  plasma  frequency 
fp  =  9  >JWt  kHz  (for  Ne  in  cm  3),  which  depends  only  on  the 
plasma  density  Ne.  High-power  radio  waves  can  transfer  sig¬ 
nificant  energy  to  the  ionospheric  plasma  when  fT  matches 
plasma  eigen-frequcncies  near  fp ,  the  upper  hybrid  fre- 

quency  Juh  =  yJ/f+jf,,  or  multiples  of  the  electron  gyro 

frequency  fce.  This  energy  can  heat  the  ionospheric  electron 
population  to  thousands  of  Kelvin  above  ambient  and  also 
accelerate  electrons  to  suprathermal  energies  in  the  range 
of  a  few  to  few  dozen  eV.  Suprathermal  electrons  excite  optical 
emissions  upon  impact  with  the  neutral  gasses  [Bernhardt 
et  al.*  1988],  and  can  actually  create  new  plasma  when  their 
energy  exceeds  the  ionization  potentials  of  the  gasses  (-12 
1 8  eV)  [Gustavsson  et  al. ,  2006]. 

[4]  Based  on  energy  budget  estimates,  creation  of  artifi¬ 
cial  ionization  by  HF  heating  has  been  predicted  to  occur 
when  power  densities  in  the  HF  beam  reach  1  GW  effective 
radiated  power  (ERP),  similar  to  that  of  the  solar  extreme 
ultraviolet  flux  creating  the  natural  ionosphere  [Carlson* 
1993].  However,  the  effects  of  ionization  production  have 
generally  been  too  small  to  be  detectable  in  the  ionosphere 
until  the  recent  upgrade  of  the  HAARP  facility  in  Gakona, 
Alaska  (62.4°  N  145°  W)  to  3  6  MW  power,  providing 
ERP  in  the  range  of  400-4000  MW  Initial  experiments  afier 
the  upgrade  showed  indications  of  enhanced  plasma  densities 
near  200  km  altitude  on  the  bottomside  of  the  ionospheric 
F  region,  which  apparently  de-focused  the  transmitter  beam 
into  a  ring  surrounding  a  bright  central  spot  visible  in  optical 
measurements  [Pedersen  et  al .,  2009].  In  this  report  we 
present  the  first  evidence  of  artificial  ionization  reaching 
the  threshold  of  self-perpetuation,  where  the  artificial  plasma 
becomes  sufficiently  dense  to  sustain  interaction  with  the 
transmitter  beam.  This  breaks  the  dependence  on  the  natural 
ionosphere  and  allows  the  artificial  plasma  to  descend  as  an 
ionization  iswave”  toward  higher  power  densities  closer  to  the 
transmitter,  forming  a  glowing  artificial  layer  near  150  km 
altitude  in  the  ionospheric  “valley”  region,  a  minimum  in 
natural  plasma  density. 

2.  Observations 

[5]  The  most  pronounced  cases  of  artificial  layer  forma¬ 
tion  in  the  valley  region  occurred  between  4-6  UT  on  March 
17,  2009,  corresponding  to  twilight  and  early  evening  hours 
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Seconds  after  05:13.00  UT  17  March  2009 

Figure  1.  (left)  Images  of  artificial  optical  emissions  as 
viewed  from  the  remote  site  at  557.7  nm  (top,  with  altitudes 
along  the  HAARP  field  line  indicated),  and  from  the 
HAARP  site  looking  up  the  field  line  with  high  resolution 
at  557.7  nm  (2nd  row)  and  427.8  nm  (3rd  row).  Average 
calibrated  intensities  at  427.8  nm  for  the  central  region  of 
the  images  are  shown  in  the  4th  row  as  a  function  of  seconds 
after  the  transmitter  turned  on  at  5:13:00  UT.  (right)  The 
557.7  nm  images  for  210s  have  been  combined  using  a  to¬ 
mographic  algorithm  to  provide  a  cross-section  of  the  opti¬ 
cal  volume  emission  rate  in  the  magnetic  meridian  plane. 
The  HAARP  magnetic  field  line  and  contours  of  nominal 
transmitter  power  in  percent  relative  to  peak  power  at 
230  km  altitude  have  been  superimposed  on  the  tomographic 
cross  section. 


in  Alaska  local  time  The  transmitter  was  pointed  at  the 
magnetic  zenith  (az  =  202°  el  =  76°)  and  operated  with 
O-mode  polarization  at  full  power  of  -440  MW  ERP  alter¬ 
nating  between  4  min  at  full  power  and  4  min  off  to  allow 
recovery  from  artificially  induced  effects.  The  transmitter 
frequency  alternated  between  3.16  and  2.85  MHz  every 
other  “on”  period  to  compare  effects  away  from  and  near 
2 fce.  After  05:05  UT  the  ionospheric  critical  frequency  fo  F2 
dropped  below  3.16  MHz  and  only  2.85  MHz  was  utilized. 
The  transmitter  beam  had  a  full  width  at  half  maximum  of 
1 8  degrees  along  the  magnetic  meridian  at  2.85  MHz.  Optical 
observations  were  carried  out  with  a  remote  wide- field  imager 
located  —  1 60  km  N  of  HAARP,  and  with  multiple  wide-  and 
narrow-field  systems  at  the  HAARP  site  observing  557.7  nm 
emissions  from  the  lS  state  of  atomic  oxygen  corresponding 
to  >4  eV  electron  energy  and  427.8  nm  Aj  emissions  indicat¬ 
ing  ionization  production  at  >1 8  eV.  Figure  1  shows  a  series 
of  optical  images  viewed  obliquely  from  the  remote  site  at 
557.7  nm  (top  row)  and  looking  up  the  magnetic  field  from 
HAARP  (557.7  and  427.8  nm  in  the  2nd  and  3rd  rows)  during 
an  artificial  layer  creation  event  from  05: 13-05:17  UT.  Emis¬ 
sion  intensities  averaged  over  the  central  part  of  the  427. 8  nm 
images  are  shown  in  the  fourth  row. 

[6]  Early  in  the  period,  the  optical  emissions  are  diffuse 
and  roughly  proportional  to  the  transmitter  beam  intensity. 
In  the  second  minute,  the  diffuse  emissions  have  organized 
themselves  into  a  sharp-edged  central  disk  with  indications 
of  a  ring  surrounding  it.  Very  bright  field-aligned  filaments 
begin  to  appear  in  the  center  of  the  disk  and  descend  in  al¬ 
titude.  In  the  3rd  minute,  the  filaments  are  very  intense  and 
the  emissions  have  descended  to  nearly  150  km  altitude.  The 


central  region  immediately  surrounding  the  filaments  also 
begins  to  grow  dark,  a  trend  which  reaches  its  extreme  in 
the  4th  minute,  where  emissions  in  the  center  of  the  beam 
have  formed  a  bright  spot  near  150  km  altitude,  leaving 
an  empty  ring  near  200  km.  The  localized  spot  and  ring 
are  clearly  shown  in  a  cross  section  through  a  tomographic 
reconstruction  generated  from  557.7  nm  images  from  the 
2  sites  (Figure  1,  right).  The  magnetic  field  line  and  nominal 
contours  of  transmitter  beam  power  in  percent  relative  to 
230  km  altitude  at  the  beam  center  have  been  superimposed 
on  the  cross  section. 

[7]  Although  the  accelerated  electrons  exciting  the  optical 
emissions  can  travel  -  10  km  or  so  along  the  magnetic  field 
from  their  source,  the  isolated  spot  near  150  km  in  the  beam 
center  indicates  that  a  strong  interaction  with  the  2.85  MHz 
transmitter  beam  is  taking  place  at  this  altitude.  Note  that  the 
natural  plasma  density  at  that  altitude  at  this  time  of  night  is 
below  2.5  x  104  cm  1  or  -200  kHz,  an  order  of  magnitude 
too  low  for  interaction. 

[8]  Ionograms  recording  the  apparent  range  of  points 
where  the  plasma  frequency  matches  the  swept  lonosonde 
probe  frequency  were  acquired  at  1  min  intervals  throughout 
the  experiment.  Ionograms  just  before  and  during  a  transmitter 
“on”  cycle  that  started  at  05:21  UT  are  displayed  in  Figure  2. 
The  “off’  ionogram  (Figure  2,  left)  shows  only  the  back¬ 
ground  F-region  ionosphere,  which  peaks  near  250  km. 
Densities  below  -200  km  virtual  height  are  below  the 
1  MHz  detection  threshold  in  the  mode  utilized.  After  sev¬ 
eral  minutes  of  heating,  however  (Figure  2,  middle),  clear 
layers  of  new  echoes  are  seen  at  -200  km  and  160  km  vir¬ 
tual  range,  closely  corresponding  to  the  observed  optical 
structures.  Although  the  HAARP  transmitter  interferes  with 
ionosonde  reception  near/7-,  X-mode  echoes  (green)  reflect 
from  the  same  density  levels  at  higher  frequencies  above 
this  gap,  and  can  be  related  to  O-mode  echoes  by  the  ex¬ 
pression/  =  y  +  \J Q  H-  (y )  =  f„  +^-.  This  indicates  that 

the  O-mode  critical  frequency  in  these  layers  reached  2.6 
and  3.0  MHz.  True  height  profiles  inverted  from  these 
traces  using  the  University  of  Massachusetts  Lowell  SAO 
Explorer  program  are  shown  in  the  right  panel,  although 
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Figure  2.  (left)  Ionogram  showing  only  background  F- 
region  echoes  while  the  transmitter  was  off.  (middle)  After 
several  minutes  of  heating,  the  interaction  region  descends 
to  lower  altitudes,  and  two  lower  layers  of  echoes  become  ap¬ 
parent  near  160  and  200  km  virtual  height,  (right)  Using  both 
the  O-mode  (red  and  pink)  and  X-mode  (green)  echoes  allows 
the  gap  resulting  from  interference  at  the  transmitter  frequency 
to  be  bridged  and  density  true  height  profiles  to  be  determined 
for  the  various  layers. 
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Figure  3.  Time-vs-altitude  plot  of  557.7  nm  optical  emissions  along  the  HAARP  field  line  with  contours  showing  the 
altitudes  where  fp  =  3.16  MHz  (blue),  fp  =  2.85  MHz  (rcd\fufl  =  2.85  MHz  (green),  and  2fce  =  2.85  MHz  (dashed  white). 
Horizontal  blips  such  as  the  one  seen  near  200  km  altitude  at  5:05  UT  are  stars  passing  through  the  view. 


the  top  sides  of  the  profiles  are  extrapolated  by  the  soft¬ 
ware  and  are  not  observable  with  the  ionosonde  technique. 
The  inescapable  conclusion  from  the  combined  optical  and 
ionosonde  measurements  is  that  artificial  ionospheric  layers 
were  produced  at  200  and  -150  km  altitude,  and  that  the 
density  in  the  lower  of  the  layers  actually  exceeded  both 
fr  and  the  background  fa  F2 . 

[9]  This  phenomenon  was  reproduced  on  multiple  nights, 
but  in  the  following  section  we  present  a  broader  view  of 
several  repetitions  of  the  experiment  on  the  night  of  March 
17.  Figure  3  shows  sequential  altitude  profiles  extracted 
from  the  oblique  optical  images  along  the  projection  of  the 
HAARP  magnetic  field  line,  overlaid  with  contours  of  plasma 
frequency  obtained  by  inversion  of  the  l -min  ionograms. 

[to]  On  this  night,  artificial  optical  emissions  at  557.7  nm 
were  first  detected  above  the  twilight  during  a  4-min  “on” 
period  at  04:49  UT  at  2.85  MHz,  but  a  gap  in  the  optical  data 
prevented  full  analysis  of  this  period.  Nevertheless,  the  iono¬ 
sonde  recorded  increased  densities  on  the  bottomside  and  in 
the  valley  region  as  evidenced  by  descending  contours  of 
fuh  on  the  bottomside  and  2.85  MHz  in  the  valley  region 
Transmissions  away  from  2 fce  at  3.16  MHz  were  attempted 
in  3  experiment  cycles  prior  to  04:45  UT,  but  produced  no 
significant  effects  and  were  discontinued  after  the  4th  cycle 
at  04:57  UT.  At  2.85  MHz  at  05:05  UT,  however,  strong 
emissions  started  out  as  a  diffuse  glow,  then  gathered  into 
the  spot-within-ring  bullseye  structure  while  gradually  des¬ 
cending  in  altitude.  Ionograms  showed  a  secondary  bottom- 
side  layer  gradually  descending,  with  the  contour  of fth 
approaching  200  km. 

[11]  The  next  “on”  period  starting  at  05: 13  UT  corresponds 
to  the  images  presented  in  Figure  l:  initial  development  is 
similar  to  the  previous  period  for  the  first  2  min.,  with  a  sec¬ 
ondary  layer  again  descending  near  200  km  at  -240  m/s.  As  the 
emissions  reached  180  km  altitude,  however,  the  557.7  nm 
intensity  suddenly  increased  by  a  factor  of  6  and  the  rate 
of  descent  rose  to  -260  m/s,  with  the  interaction  region  des¬ 


cending  to  -150  km  altitude  over  the  next  90  seconds  or  so. 
Here  the  ionograms  indicated  the  presence  of  plasma  density 
up  to  3.1  MHz,  well  in  excess  0 f/r.  During  the  4th  minute, 
the  emissions  appear  to  quench  themselves  somewhat  and 
retreat  in  altitude  before  the  end  of  the  transmission.  During 
a  continuous  “on”  period  beginning  at  05:2 1UT,  the  same 
pattern  is  realized,  with  the  emissions  brightening  at  about 
180  km  and  rapidly  descending  to  -150  km.  In  this  case 
the  emissions  appear  to  quench  themselves  several  times, 
initiating  the  descending  process  over  again  from  higher  alti¬ 
tudes,  although  some  of  the  apparent  vertical  motion  in  the 
optical  data  toward  the  end  of  the  period  represents  horizontal 
displacement.  The  ionosonde  data  indicate  that  artificially 
produced  plasma  density  in  excess  of  the  2,85  MHz  transmit¬ 
ter  frequency  was  present  most  of  the  time  optical  emissions 
were  seen  near  150  km  altitude,  as  shown  by  the  red  contours 
between  150  and  160  km  altitude.  Once  the  background 
ionosphere  had  decayed  below  fuh  at  about  05:40  UT, 
production  of  both  optical  emissions  and  artificial  iono¬ 
spheric  layers  ended.  Just  prior  to  the  end  of  the  optical  emis¬ 
sions,  the  authors  present  at  HAARP  during  the  experiment 
went  outside  and  were  able  to  observe  the  artificial  optical 
emissions  with  the  naked  eye,  an  indicator  of  the  extraordi¬ 
nary  intensity  of  the  emissions  produced  [cf.  Pedersen  and 
Gerken ,  2005], 


3.  Discussion 

[12]  We  interpret  these  observations  as  follows:  near  the 
double  resonance  where  fT  simultaneously  matches  both 
fuh  and  2 \fce,  electrons  are  efficiently  accelerated  [e  g.,  Kosch 
et  al ,  2007]  and  ionize  the  neutral  gas,  adding  to  the  ambi¬ 
ent  plasma  density.  This  pushes  the  density  contours  down 
over  a  few  minutes  to  form  a  descending  artificial  ionospheric 
layer  near  200  km,  sufficiently  dense  to  sustain  interaction  at 
fuh  but  not  fp.  This  part  of  the  process  was  seen  in  each  of  the 
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2.85  MHz  cases  on  March  1 7  as  well  as  earlier  experiments, 
and  appears  in  ionograms  as  a  secondary  bottomside  layer. 

[n]  As  the  interaction  altitude  descends,  the  power  densi¬ 
ty  in  the  transmitter  beam  goes  up  as  already  reaching 
132%  at  200  km  compared  to  230  km.  This  process  con¬ 
tinues  until  about  180  km  (163%  power),  when  the  plasma 
frequency  in  the  descending  layer  reaches  fr-  This  would 
correspond  to  formation  of  the  bright  field-aligned  structures 
in  the  optical  images.  Stimulation  of  the  plasma  resonance 
now  becomes  possible,  providing  an  additional  potential 
mechanism  for  accelerating  even  more  electrons.  Propaga¬ 
tion  of  transmitter  power  to  altitudes  above  the  artificial  layer 
is  now  blocked  in  the  center  of  the  beam,  with  any  HF  power 
not  absorbed  being  reflected  back  down  and  adding  to  the 
field  strength  below  the  reflection  point.  With  this  new  elec¬ 
tron  source,  the  localized  artificial  plasma  now  becomes  com¬ 
pletely  self-sustaining,  and  rapidly  propagates  downward  as 
an  ionization  wave  front  toward  the  transmitter,  appearing  as 
a  layer  of  echoes  in  the  ionograms  near  150  km.  For  a  real¬ 
istic  density  gradient  of  105  cm-3  over  10  km,  the  observed 
propagation  implies  an  ionization  production  rate  in  excess 
of  losses  by  about  2.6  *  103  cm-3  s_l.  The  descent  of  the 
artificial  ionization  stalls  at  ^  150  km  altitude,  near  the  tran¬ 
sition  from  atomic  oxygen  to  short-lived  molecular  ions, 
which  decay  at  rates  exceeding  even  the  enhanced  ionization 
production  rate,  now  driven  by  235%  nominal  transmitter 
power. 

[14]  The  complete  quenching  of  the  artificial  ionization 
upon  reaching  —  1 50  km  altitude  is  more  difficult  to  explain, 
but  as  seen  in  the  427.8  nm  intensity  (Figure  1 ,  bottom),  the 
overall  ionization  rate  averaged  over  the  central  region  of 
the  beam  appears  to  go  down  steadily  as  the  central  region 
descends  in  altitude,  and  the  matching  frequency  for  2 f  .€  in¬ 
creases  to  -2.95  MHz  at  150  km  altitude.  Excitation  near 
the  2fce  resonance  has  a  significant  frequency  dependence 
[Mishin  et  aiy  2005;  Kosch  et  ai,  2007]  This  likely  reduces 
the  efficiency  of  electron  acceleration  even  at  higher  power 
levels  and  contributes  to  the  artificial  plasma  density  produc¬ 
tion  balance  changing  over  from  a  net  surplus  evidenced  by 
downward  motion  of  the  layer  to  a  net  loss,  which  rapidly 
leaves  it  unable  to  interact  with  the  transmitter  beam  and 
maintain  itself.  This  can  be  tested  in  future  experiments  by 
increasing  the  transmitter  frequency  as  the  layer  drops  to 
maintain  fT  near  2 fce.  In  any  case,  as  seen  in  the  period  after 
5:21  UT  when  the  iransmitter  remained  on  continuously,  the 
interaction  quenches  itself  after  a  few  minutes  at  1 50  km  and 
starts  over  again  near  fufl  in  the  background  ionosphere. 
When/,/,  is  no  longer  available  due  to  decay  of  the  back¬ 
ground,  the  process  is  no  longer  able  to  restart  itself 

[15]  We  note  that  while  we  have  discussed  evolution  of 
the  artificial  ionization  in  terms  of  altitude,  there  is  a  parallel 
spatial  evolution  occurring  simultaneously,  with  the  bright 
filaments  appearing  away  from  the  beam  center  and  gener¬ 
ally  moving  inward  before  extinguishing  themselves.  Also, 
although  we  have  referred  to  the  regions  of  artificially  en¬ 
hanced  plasma  density  as  “layers,”  based  on  their  appear¬ 
ance  in  the  lonosonde  data,  the  optical  images  clearly 
show  a  high  degree  of  spatial  structure,  with  the  areas  of 


greatest  density  enhancement  likely  located  inside  the  very 
bright  field-aligned  filaments  a  few  km  in  diameter. 

[16]  While  the  exact  mechanisms  producing  the  artificial 
ionization  layers  in  these  experiments  remain  to  be  deter¬ 
mined  and  the  patchy  and  unstable  nature  of  the  enhance¬ 
ments  pose  difficulties  to  be  overcome  before  practical 
applications  become  feasible,  these  experiments  have  dem¬ 
onstrated  the  first  production  of  self-sustaining  artificial 
ionospheric  plasma  and  represent  a  major  breakthrough 
laying  the  groundwork  for  a  variety  of  controlled  near-space 
radio  applications  based  on  a  new  technology  of  artificial 
ionospheric  plasma  production. 
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